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We present results for one-loop matching factors of four-fermion operators composed of HYP- 
smeared staggered fermions. We generalize previous calculations by using the tree-level improved 
Symanzik gauge action. These results are needed for our companion numerical calculation of Bk 
and related matrix elements. We find that the impact on one-loop matching factors of using the 
improved gluon action is much smaller than that from the use of either HYP smearing or mean- 
field improvement. The one-loop coefficients for mean-field improved, HYP-smeared operators with 
the Symanzik gauge action have a maximum magnitude of O(l) x a s , indicating that perturbation 
theory is reasonably convergent. 
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I. INTRODUCTION 

Numerical simulations of lattice QCD are now able to 
calculate a range of phenomenologically interesting non- 
perturbative quantities with high precision. Of particular 
interest are hadronic matrix elements of operators that 
appear in the electroweak Hamiltonian, or in extensions 
of the standard model. For such quantities it is necessary 
(in order to make use of Wilson coefficients calculated in 
continuum perturbation theory) to determine the match- 
ing factors which relate operators regularized on the lat- 
tice with those regularized in the continuum. For the 
latter one typically uses naive dimensional regularization 
(NDR) with MS subtraction. 

In this paper we calculate matching factors (which are, 
in general, matrices) for four-fermion operators composed 
of light staggered quarks. These arise, for example, in the 

calculation of the K — K mixing parameter Bk- Their 
flavor structure forbids mixing with lower-dimensional 
operators, so the matching is only between operators of 
dimension 6. In the electroweak theory, the operator that 
arises has a "left-left" structure, due to the left-handed 
coupling of the W-bosons. In extensions of the standard 
model, however, AS" = 2 operators can arise with other 
Dirac structures. For this reason we calculate matching 
factors for all possible Dirac structures. 

In recent years, it has become increasingly common to 
determine matching factors non-perturbatively, either us- 
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ing the Rome-Southampton non-perturbative renormal- 
ization method [![, or using approaches based on the 
Schrodinger functional Q. These methods replace hard- 
to-estimate truncation errors by controllable statistical 
and systematic errors. We are implementing such calcu- 
lations for improved staggered fermions, but have so far 
only obtained results for bilinear operators Q. We ex- 
pect that the implementation for four-fermion operators, 
which involves mixing with a long list of lattice operators, 
will be more challenging. The use of one-loop matching 
is a useful intermediate step, and, as we will describe, 
the necessary calculations are relatively simple general- 
izations of previous work. We also note that, since our 
present numerical calculations involve very small lattices 
spacings (a 0.045 fm), the truncation errors are quite 
small, since they are proportional to a 2 s with a s evaluated 
at a scale ss 1 a [4j. 

Our companion numerical calculations use valence 
staggered fermions which have been improved by the use 
of HYP-smeared links (links replaced with hypercubic 
blocked links Q). The ensembles are those generated 
by the MILC collaboration [6], in which the gauge ac- 
tion is Symanzik-improved, and the quark action is the 
asqtad staggered action. Previous calculations have ob- 
tained the matching factors for four-fermion operators 
composed of HYP-smeared staggered fermions Q, but 
only using the Wilson gauge action. Here we generalize 
these results to the case of an improved gauge action. 
This extends our earlier work in which we calculated 
matching factors for bilinear operators using improved 
gauge actions and HYP-smeared staggered fermions [8[. 

At first sight, the generalization from the Wilson gauge 
action (for which, in Feynman gauge, the gluon propaga- 
tor is diagonal in Euclidean indices) to an improved gauge 
action (in which the propagator is not diagonal) appears 
non-trivial. In particular, one-loop calculations using 
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unimproved staggered fermions 0, Eo| were simplified us- 
ing the diagonal nature of the gluon propagator. The 
inclusion of smeared links, however, leads to the natural 
introduction of a "composite gluon propagator" which 
represents both the smearing and the gluon propagator 
itself. This propagator contains non-zero off-diagonal el- 
ements, and so calculations of one-loop matching factors 
for HYP-smeared staggered fermions must already deal 
with the presence of such elements @, [H| ■ This means 
that the generalization to an improved gluon propaga- 
tor requires no change to the analytic expressions — all 
one needs to change is the composite gluon propagator 
before numerical evaluation of the loop integral. As we 
discuss here, this simplification holds not only for bilinear 
operators [8j, but also for four-fermion operators. 

The paper is organized as follows. In Sec. HU we re- 
call our notation and conventions for actions and opera- 
tors. In Sec. Mil we present the Feynman diagrams and 
describe their evaluation. Because we are building on 
the work of Refs. 0] and @, we provide only a minimal 
discussion of technical details. In Sec. IIV1 we present 
our numerical results for matching factors, providing the 
complete matching matrix for the operators relevant to 
Bk 1 and a partial matrix (the part that will likely be 
used in practice) for other four-fermion operators. We 
conclude briefly in Sec. fVl 
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Here, h llu {k) is the smearing kernel, which depends on 
the smearing parameters and the details of the HYP con- 
struction. It contains non-zero off-diagonal components 
because a smeared link in one direction contains contri- 
butions from thin links in all four directions. It turns 
out that we need only the linear term in Eq. §5§ in a 
one-loop calculation. The contribution of the quadratic 
term (which gives rise to tadpole diagrams) turns out to 
vanish due to the projection back into the SU(3) group 
that is part of the definition of HYP-smearing [9|, 53 ■ 
The smearing kernel h^ u can be conveniently decom- 
posed into diagonal and off-diagonal parts: 

h(u,(k) = S^D^k) + (1 - S^s^Gu^k) , (6) 
with s M = sin(fc M /2), and 
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II. ACTIONS AND OPERATORS 

The HYP-smeared staggered action has the same form 
as the unimproved staggered fermion action, 
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where rj^ri) = (— l)" lH ^'V-i ; anc ] the covariant differ- 
ence operator is 

V**(n) = hy^xin + A) - V}(n - fa) X (n - fa)} • (2) 

Here and in the following we set the lattice spacing a to 
unity, except where confusion could arise. HYP improve- 
ment consists of using HYP-smeared links, Va, instead of 
the original "thin" links, Up. We set the HYP-smearing 
parameters to the values that remove the tree-level cou- 
pling of quarks to gluons having one or more components 
of momenta equal to 7r. These values are ct\ = 0.875, 
a 2 = 4/7 and — 0.25, in the notation of Ref. @. 
These are the values used in our numerical simulations. 

After gauge- fixing, we expand both the thin and 
smeared links in the usual way, 



U^(n) = expligoA^n + (1/2)] 
V^(n) = exp[ig -Bp(n + fa/ 2)] 



(3) 
(4) 



where go is the bare gauge coupling. The relation be- 
tween the fluctuations of smeared and thin links can be 



Here /i, is, p, and a all differ from each other. By contrast, 
the smearing kernel for an action containing the original 
thin links is simply = 

We use the tree-level Symanzik-improved gluon action 
El El; 
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where "pi" and "rt" represent plaquette and rectangle, 
respectively. In fact, the MILC collaboration use the 
(partial) one-loop Symanzik-improved action determined 
in Refs. [HI, [T3|- However, the one- loop contributions 
to this action contribute to matching factors of valence 
fermionic operators only at two-loop level, so the consis- 
tent choice for our one-loop calculation is the tree-level 
action ©. For purposes of comparison, we also use the 
Wilson gluon action, which is obtained from Eq. © by 
dropping the rectangle term and setting the coefficient of 
the plaquette to unity instead of 5/3. 

Since we use MILC asqtad ensembles in our numeri- 
cal studies, the sea quarks are asqtad staggered fermions 
rather than HYP-smeared. We do not display the sea- 
quark action, however, since sea-quarks only enter at 
two-loop order in the matching of valence fermionic op- 
erators. Our one- loop matching factors are thus valid for 
any choice of sea quarks. 

We now turn to the definitions of our lattice four- 
fermion operators, which are the same as those used in 
Ref. [7J] . Our construction follows the hypercube method 
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of Ref. [18j. The operators come in two classes, differing 
in the contractions of their color indices. First we have 
one color-trace operators, labeled with a subscript /: 

[S X F][S' x F'}j(y) = 

& E [x^^y + A){ ls ®i F ) ABX f\2y + B)] 

A,B,A',B' 

x [£\2y + A')T^®^) A , B , X ^\2y + B')] 

x V ab '(2y + A,2y + B')V a ' b (2y + A', 2y + B) . (10) 

Here, y € Z 4 is the coordinate of 2 4 hypercubes. Hyper- 
cube vector^] S and S' denote the spins of the component 
bilinears, while F and F' denote the tastes. Indices a, 
6, a', and b' denote colors, while superscripts (i) for i = 
1, 2, 3, 4 label different flavors (not tastes). Using four dif- 
ferent flavors forbids penguin diagrams, which would lead 
to mixing with lower-dimension operators!^ Two "fat" 
Wilson lines V ab ' (2y+A, 2y+B') and V a ' b (2y+A', 2y+B) 
ensure the gauge invariance of the four-fcrmion operators. 
A fat Wilson line V ab ' (2y + A, 2y + B'), for example, is 
constructed by averaging over all the shortest paths con- 
necting 2y + A and 2y + B', with each path formed by 
products of HYP-smeared links V^. When we use the 
unimproved staggered action the Wilson lines are com- 
posed of unsmeared thin links, U^. 

The second class are the two color-trace operators, for 
which we use the subscript II: 

[S x F][S' x F'] u (y) = 

1 J2 [x i a 1) ^y + A)J^Ii) AB X ( " ) (2y + B)} 

A,B,A',B' 

x [xf{2y + A'){ ls ,®Z F ,) A , B , X %\2y + B')] 

xV ab (2y + A,2y + B)V a ' b ' (2y + A', 2y + B') . (11) 

These operators differ from those with one color-trace 
only by the choice of fat Wilson lines — here they con- 
nect within each bilinear, whereas for the one color-trace 
operators they connect between bilinears. 

We also consider mean-field improvement of the stag- 
gered action and operators following Refs. [1, [l9l - l2lj . 
This is also referred to as tadpole improvement. Mean- 
field improvement is achieved by rescaling the staggered 
fields and the links. In the case of HYP-smeared stag- 
gered fermions the rescaling is 
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with V p \ the plaquette composed of HYP-smeared links. 
One then constructs the operators described above out 
of ip, ip, and Vfj,. The resulting operators are expected to 
be closer to their continuum counterparts because the 
rescaled links fluctuate around unity. Mean-field im- 
provement can be implemented in simulations after the 
data has been collected, as long as the contributions to 
the four-fcrmion operators having different numbers of 
links are stored separately. 



III. FEYNMAN DIAGRAMS AND THEIR 
EVALUATION 

Feynman rules for the gauge and staggered- fermion ac- 
tions, and for insertions of the four-fermion operators, 
can be found in literature and we do not reproduce them 
here. The rules for unimproved staggered fermions are 
given in Refs. P, ll9L 22], and the generalization to HYP- 
smeared staggered fermions can be found in Ref. 0, EH ■ 
The gluon propa gat or for the Symanzik action was de- 
termined in Ref. [14fl; we use the simpler form presented 
in our earlier work 12]. 

We show the Feynman diagrams contributing to one- 
loop matching factors to the two types of four-fermion 
operators in Figs. [T] and [2] Analytic formulae for these 
diagrams for HYP-smeared staggered fermions with the 
Wilson gluon action are given in Ref. @. We do not 
repeat these results here, since it turns out, as already 
mentioned in the Introduction, that the generalization to 
the improved gluon action is relatively simple. Instead 
we explain the recipe by which the results of Ref. Q can 
be generalized. 

The key point is that, since all gauge links are HYP- 
smeared (whether in the action or the operators), the 
gluon propagator always comes with a smearing kernel 
on each end. Thus what appears is the composite gluon 
propagator (called the "smeared-smeared propagator" in 
Ref. |l3): 



T MV (k) = E h ti.a{k)h v p{k)'D a i } {k) . 



(16) 



a/3 



Here /x and v are the directions of the initial and final 
smeared gauge links, h is given in Eq. (j6]), and T> a p(k) is 
the gluon propagator in Feynman gauge. Even with the 
Wilson gauge action, where T> is diagonal, the fact that 
h has non-vanishing off-diagonal elements implies that 
T does too. Thus the generalization to the Symanzik 
gauge action, for which T> itself has non-vanishing off- 
diagonal elements, does not introduce any fundamentally 
new types of contribution to T . Of course, the expression 
for T is much more involved, but this does not present 



1 These are vectors whose entries are or 1. 

2 The relation of these four-flavor operators to the AS ■■ 
uum operators is discussed below. 



3 To be precise, we use the formulae of Appendix A of Ref. [T^| 
2 contin- with u = 1, c = —1/12 and c' = 0. The result for the Wilson 

gauge action is obtained by further setting c = 0. 
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FIG. 1. Feynman diagrams contributing to the one-loop ma- 
trix elements of one color-trace operators. We show only one 
diagram of each type. Hypercube vectors (C, D, C' , and D 1 ) 
multiplied by n denote external quark momenta, a, b, a' , and 
b' are color indices. Dashed lines indicate the Wilson lines 
which make the four-fermion operator gauge invariant. Boxes 
indicate the hypercube bilinears of which the four-fermion op- 
erator is composed. 



FIG. 2. Feynman diagrams contributing to the one-loop ma- 
trix elements of two color-trace operators. Notation is as in 

Fig.ru 
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problems since the expression is evaluated numerically 
when doing the loop integral. This situation is in con- 
trast to what happens if the action and operators are 
composed of thin links, for then T is diagonal, which 
greatly simplifies the resulting expressions. 

In order to simplify the expressions for Feynman dia- 
grams, Ref. used the following properties of T^ v : it 
is symmetric, and its off-diagonal elements are propor- 
tional to s M s„ multiplied by a function even in each of 
the components of For the Wilson gauge action, with 
diagonal T>, these properties follow from the fact that 
has the same properties. For the Symanzik gauge action, 
it turns out that T> also has these properties, from which 
it is simple to show that T does too. Thus the simplifi- 
cations used in Ref. apply for both gauge actions. 

We now describe how the analytic formulae of Ref. [7[ 
must be changed when using the Symanzik action@ Two 
independent approaches to the matching calculation were 
used in that work. In the first, explicit expressions 
were given for all diagrams (Appendices A, B and C of 
Ref. 0). To obtain the expressions for the Symanzik 
gauge action one must make the following replacement: 



A 



a/3 



Imp 



(17) 



In the second method (Appendix D of Ref. [7] ) , maximal 
use was made of the matching calculation for bilinear 
operators. For this part of the calculation, one can simply 
use our results for matching factors of bilinear operators 
with the Symanzik gauge action [8j. For two classes of 
diagrams [those of Figs. [5Jf) and (h)], bilinear results 
are not sufficient, and for these Ref. gives explicit 
expressions. These are written in terms of the diagonal 
and off-diagonal parts of T, and in particular in terms of 
and O m „ defined through 



5^ V P^ + (1 - S^As^s^Of, 



4E 



(18) 



(where repeated indices are not summed). Here the 
recipe is to construct T^ v , Eq. (TT()|) . using the Symanzik 
gluon propagator, use this in Eq. (fT8"| to obtain new ex- 
pressions for and , and use the latter in the results 
of Ref. 0. 

As in Ref. , we evaluate matching factors using both 
methods described above and find agreement. This is a 
non-trivial check on the numerical implementation of the 
analytic expressions. We have also checked the relations 
which follow from Fierz identities and from the U(l) e 
symmetry of staggered fermions. 



IV. MATCHING FACTORS 

We calculate the matching factors in the usual way 
by evaluating the qqqq matrix elements of the operators 
both on the lattice and in the continuum, and projecting 
onto the different color and spin-taste contributions. We 
do so at one-loop order, requiring the evaluation of the 
diagrams of Figs. Q] and [2] on the lattice. On the contin- 
uum side, only diagrams of types (a) and (g) contribute, 
since the continuum four-fermion operators do not con- 
tain gauge fields. In the continuum calculation one must 
choose the continuation to 4 + e dimensions of the op- 
erators and the projectors onto different spin structures. 
We follow the conventions described in Refs. and [23| . 

The matching formula between continuum and lattice- 
regularized operators then takes the general form 



O 



Cont 



( / i) = ^Z ii ( M) a)0 J Lat (l/a), 



(19) 



with /i the continuum regularization scale, and the lattice 
spacing now made explicit. At one-loop order, and with 
a suitable choice of lattice operators, the matching factor 
has the form 



(4tt) 2 



0(a). (20) 



where 7^ and cy are, respectively, the one-loop anoma- 
lous dimension matrix and the finite coefficients. The 
latter are given by the difference of finite terms in the 
continuum and lattice one-loop calculations, 



C 



Cont 



c 



Lai 



(21) 



The general expressions for 7,-,- and Gy ont are given in 
Ref. and we do not reproduce them here@ We only 
note that the mixing structure in the continuum is much 
simpler than that on the lattice because taste is con- 
served. 

Mean-field improvement of the action and operators 
leads to a change in C^ at and thus in the finite part of 
the matching factors: 



MF 



Cij — CfImfTv, 



where Cp = 4/3, and 



IMF 



(4tt)< 



L(2^( (s " 2)271l - S " lS " 2Tl2 



(22) 



(23) 



4 As discussed in Ref. Q , the simple recipe described here does not 
work if one uses the asqtad action because not every diagram can 
be expressed in terms of the composite gluon propagator (due to 
the presence of the Naik term). Some diagrams would need to 
be calculated anew. 



5 In Ref. 0, a more elaborate notation is used in which 7 and C 
become matrices in "color-trace" space. We do not use this nota- 
tion here. We take this opportunity to correct two typographical 
errors in Table XIV of Ref. Q: the entries in the 7^ column 
which are —4 and 4/3 should be changed to +4 and —4/3, re- 
spectively. 
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with Tin, the composite gluon propagator defined in 
Eq. (|16[) . General results for the mean- field-improvement 



coefficients, T^-, can be found in Ref. 20], and are quoted 
below for the operators considered here. 

For each of the indices i and j, there are 16 4 choices 
of the S, F, S', F'. Although lattice symmetries reduce 
the number of independent entries, remains a large 
matrix. We have obtained expressions for all its entries, 
but present here only the most interesting subset. 



A. Matching Factors for Bk 

The continuum AS = 2 four-fermion operator whose 
matrix element enters into the kaon mixing parameter 
B K is 

C bT = " 75K] \s h l»{l - To)d»] • (24) 

In order to calculate (A' |C§° nt |ifo) using staggered 
fermions, one must first relate it, in the continuum, to 
a matrix element in an augmented theory in which there 
are four tastes for each continuum flavor. In fact, as ex- 
plained in Ref. [24| , one also needs to choose the quarks in 
each bilinear to have different flavors [as has been done 
in the lattice operators defined in Eqs. (fTtH and (fTTj) ]. 
This is necessary so that the Wick contractions in the 
original and augmented continuum theories agree. Thus 
one ends up with an eightfold increase in the number 
of valence flavors. To maintain the equivalence of the 
sea-quark sectors (and thus the dynamics) of these two 
theories one must take the 8th root of the fermion deter- 
minant. This rooting is not controversial in the formal 
continuum limit, and the equality of the corresponding 
matrix elements in the two continuum theories is valid 
non-perturbatively. We do note, however, that the need 
for rooting in the augmented theory implies that this the- 
ory is partially quenched, as has been stressed in Ref. [1^] . 

The end result, in the "two spin-trace" formulation of 
Refs. [13, l2a |. is that the relevant operator in the aug- 
mented continuum theory is (keeping only the positive 
parity part) 



OCont' /T-jCont/ 
B K — U V1 
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(26) 
(27) 
(28) 
(29) 



Here S, D, S' and D' are Dirac fields having an im- 
plicit taste index running over four values. This index 
is contracted with the taste matrix £5. The overall nor- 
malization of this operator is unimportant as it cancels 
in the ratio which defines Bk- Note also that the choice 
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of taste matrix is arbitrary in the continuum theory — we 
use £5 since that is what is used in lattice calculations. 

The augmented continuum theory has been chosen to 
be the continuum limit of the lattice staggered theory0 
In particular, at tree level, the operators 0<? ont match 
onto lattice operators [defined in Eqs. (fTU|) and (fTTj) ] as 
follows: 



(30) 
(31) 
(32) 



Thus the tree-level matching relation for the Bk operator 
is 

^Cont' = ^Lat + ^Lat + pLat + ^Lat + Q ^2j + 

(34) 

At one-loop order, many lattice operators contribute 
to the matching formula. It is convenient to divide them 
into the two classes: (A) the four operators which arise 
at tree-level, which have the £5 taste matrices in both bi- 
linears, and (B) the remaining operators, which all turn 
out to have taste matrices other than £5 in the bilin- 
ears. This division is useful for two reasons. First, in 
present numerical calculations only operators from class 
(A) are kept, so these are the matching coefficients that 
are neededQ Second, these are the only operators for 
which anomalous-dimension matrix elements and finite 
continuum coefficients are non-zero. Thus we write the 
one-loop matching formula as 



-jCont 
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ie(A) 
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Lai 
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_/Lat 



(36) 



where the subscripts to the sums indicate that i runs 
over the four operators in class (A) while j runs over 
all operators in class (B). We have put in the values of 
the anomalous dimensions and finite coefficients. The 
constants at are obtained by summing elements of the 



(25) matrices Cjf* introduced above. 



Numerical values for the are given in table |U 

We compare the naive staggered action (with operators 
having thin links) to the HYP-smeared staggered ac- 
tion (with operators having smeared links). In the for- 
mer case, we implement mean-field improvement (since, 



6 Here we assume that rooting introduces no problems with the 
continuum limit, following the discussion in Rets. 25, 27-29]. 

7 The rationale for this is that we use external kaons with taste 
£5. As shown in Ref. [30|. however, leaving out the operators 
with other tastes leads to an error of 0(a s m? K /Aj), which is 
of next-to-leading order in staggered chiral perturbation theory. 
This error must be accounted for when fitting. 
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TABLE I. Results for d^ at for various choices of gauge and fermion action and of the four-fermion operators. 





(a) 


(b) 


W 


(d) 


00 


(/) 


Gluon action 


Wilson 


Sym 


Wilson 


Sym 


Wilson 


Sym 


Quark action 


Naive 


Naive 


HYP 


HYP 


HYP 


HYP 


Mean-field imp. 


Y 


Y 


N 


N 


Y 


Y 


7Lat 
"VI 

7Lat 
"V2 

7Lat 
"Al 

7Lat 
a A2 


-2.349(1) 
-12.915(2) 
-2.951(1) 
-3.725(1) 


-2.487(1) 
-11.537(2) 
-3.077(1) 
-2.895(1) 


-4.984(1) 
-11.108(2) 
-5.496(1) 
1.012(1) 


-3.649(1) 
-8.584(2) 
-4.119(1) 
1.087(1) 


-2.174(1) 
-5.487(2) 
-2.686(1) 
1.012(1) 


-1.722(1) 
-4.729(2) 
-2.192(1) 
1.087(1) 


Range 


10.57 


9.04 


12.18 


9.67 


6.50 


5.82 



in general, perturbation theory is very poorly conver- 
gent without this improvement), while for the HYP- 
smeared case we show results both with and without 
mean-field improvement. For each choice of fermion, we 
compare the results obtained using the Wilson and tree- 
level Symanzik gauge actions, with the former being from 
Ref. [7[. We see that improving the gauge action has a 
small effect, which, in most cases, reduces the size of 
the coefficients. A much more significant reduction is 
obtained by HYP smearing, as can be seen by compar- 
ing, for example, columns (b) and (f). (This is the ap- 
propriate comparison because both columns show results 
with mean-field improvement implemented.) The results 
needed for our companion numerical calculation [3l| are 
those of column (f). 

As can be seen from Eq. (j3"6")l . the full one- loop correc- 
tion includes the anomalous dimension and finite contin- 
uum contributions as well as d^ at . Thus the total size of 
the one loop correction depends on the renormalization 
scale and lattice spacing through the combination \ia. A 
better measure of the size of the correction is the range 
of the coefficients d^^, from which anomalous dimen- 
sion and continuum contributions cancel. The ranges are 
given in table U and show a small reduction with the use 
of the improved gauge action. 

To give a sense of the numerical size of the matching 
coefficients themselves, we show in table [TT1 results for the 
zi-4 for the "ultrafine" MILC lattices (a w 0.045 fm), 
setting fi — 1/a ( "horizontal matching" ) , and using a s = 
.9 2 /(4tt) = 0.2096 (the value in the MS scheme at fi = 
1/a). For the actions we use in practice [column (f)] the 
one-loop corrections range between +2% and —8%. 



improving the gauge action leads, as above, to a small 
reduction in the magnitude of all the matching coeffi- 
cients. Note that, since these mixings are pure lattice 
artifacts, with no anomalous dimension or other contin- 
uum contributions, reducing the size of the coefficients is 
an unambiguous improvement. 

As noted above, in present numerical calculations the 
mixing with these operators is not being incorporated in 
the lattice operators, but rather is a source of systematic 
error that must be estimated by fitting. An alternative 
approach would be to include the dominant operators 
from the Tables in the one- loop matched operator. As 
one can see, there are relatively few operators having 
0(1) coefficients: 

1. [Sxvjsxyj, 

2. [P x Vp)[P x and [P x V^P x V„], 

3. Pjji/ x VjpV x ly, [T^ x V P ][T^ x V p \ 
[T„ v x V„][T„ V x V v ] and [T^ x V p ][T^ x V v ], 

The remainder of the coefficients are an order of mag- 
nitude or more smaller (i.e. K Lat | < 0.2). Since these 
coefficients are multiplied by g 2 /(47r) 2 sa 0.017 - 0.025 
for a w 0.045 — 0.12 fm, we expect the contributions to 
Bk from the remaining operators to be very small. 



Matching Factors for other four-fermion 
operators 



TABLE II. Values of Zj for the MILC ultrafine ensembles. 
The notation for actions is as in table HI 





(a) 


(b) 


(c) 


(d) 


(e) 


(/) 


ZV1 


0.978 


0.980 


1.022 


1.000 


0.975 


0.968 


ZV2 


1.154 


1.131 


1.124 


1.082 


1.030 


1.018 


ZA\ 


0.988 


0.990 


1.031 


1.008 


0.984 


0.975 


ZA2 


1.001 


0.987 


0.922 


0.921 


0.922 


0.921 



To give a complete view of the one-loop matching, we 
present, in Tables ITTI1 and ITVl the coefficients d^ at for all 
other operators which appear at this order. We see that 



Models of new physics can lead, after integrating out 
heavy particles, to AS — 2 operators with different Dirac 
structure from that in To constrain these mod- 

els one needs to know the matrix elements of these new 
operators. A standard basis is 32] 



(37) 
(38) 
(39) 
(40) 





[s a (l 


-75Kp fc (l 


-75)0* 


^Cont = 


[s a (l 




- 7 5 k; 


^Cont = 


[s a (l 


-7sKp fc (l 


+ 75)0* 


^Cont 


[s a (l 


-7 5 )d b ][s b (l 


+ 75K 
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TABLE III. Matching coefficients dj at [defined in Eq. 
for the operator required for calculating Bk, and for oper- 
ators j which have different taste than the continuum oper- 
ator (|25p . Lattice operators and fermion action are HYP- 
smeared and the gauge action is either Wilson — column (c) — 
or Symanzik — column (d). Results in column (c) are obtained 
from Tables I- IV of Ref. 0. The coefficients Tj give the im- 
pact of mean-field improvement: dj at — > d^ at + CfImfTj, 
with Imf = 1.053786 for the Wilson gauge action and 
Imf = 0.722795 for the Symanzik gauge action. Greek in- 
dices are implicitly summed, with the condition that they 
are unequal, and the further constraint that for the oper- 
ator [Vp x T„ p ][Vp x T„ p ], v < p, while for the operators 
x V P ][T^ x V p ] and [T M „ x V P ][T P „ x V v ], p. < v. Re- 
sults are accurate to at least ±2 in the last digit quoted. 





color trace 




(d) 


T 


[U A Vpj[0 A V pj 


j 


—3.450 


—2.805 


I 


\S x V, 1 \S x V ] 

yu a v pj yu A V pj 


n 


—0.263 


—0.249 





S X V S" X V 1 


I 


0.118 


0.108 


o 


,S X V S X V 
[iJJ A V pj [*J A VivJ 


H 


—0.104 


—0.097 


o 


9 x A US' x /I 1 


I 


0.043 


0.028 


o 


S X A S X ,4 


H 


—0.052 


—0.035 


o 


S X A S X /I 1 


I 


—0.015 


—0.010 


o 


Lf x A S x ,4 


H 


0.002 


0.002 


o 


rv x SUV x SI 

yv p a j yyv p a 


I 


—0.044 


—0.029 


o 


rv x s irv x s 

[V (1 A [ Vp A J J 


H 


—0.008 


—0.005 


o 


V x T HV x T 1 
[ Vp a ± hi/ yyv p a ± pyj 


I 


—0.124 


—0.086 


o 


V x T V x T 

[ Vp A J- hi/ \yv fj, A J pi/J 


H 


—0.114 


—0.084 


o 


v, x t„„ nv, x T„J 


I 


0.029 


0.023 


o 


[V/i x t m „][v m x T MP ] 


II 


-0.023 


-0.019 





[V M x Tp„][Vp x T„ p ] 


I 


0.016 


0.014 





[v„ x Tp„][Vp x r„ p ] 


II 


0.002 


0.002 





[v„ x T,, p ][Vp x r^j 


I 


0.016 


0.014 





[v„ x T,, p ][Vp x r^j 


II 


0.002 


0.002 





[v^ x r„ p ][Vp x r„ p ] 


I 


-0.118 


-0.091 





[v M x t vp \[v^ x t„ p ] 


II 


-0.037 


-0.030 





[V M x T,, p ][Vp x T^] 


I 


0.027 


0.022 





[v„ x T,, P ][V M x n,,] 


II 


-0.020 


-0.016 





[T M „ x Vp.][T M „ x V M ] 


I 


2.071 


1.547 


-1 


[T„„ x V M ][T M „ x Vj»] 


II 


-0.538 


-0.485 





[T„„ x Vfi] [T M „ x V„] 


I 


-0.410 


-0.383 





[Tp„ x Vp.][T M „ x V„] 


II 


0.452 


0.417 





[T M „ x Vp.][T M „ x Vp] 


I 


0.129 


0.120 





[T M „ x V M ][T M „ x V p ] 


II 


0.126 


0.118 





[T M „ x V P ][T M „ x 1/ M ] 


I 


0.129 


0.120 





[T M „ x V p ][Tp„ x Vp] 


II 


0.126 


0.118 





[Tp„ x V p ][Tp„ x V p ] 


I 


-2.930 


-2.331 


1 


[Tp„ x V p ][Tp„ x V p ] 


II 


-0.346 


-0.316 





[Tp„ x Vp][Tp„ x V,] 


I 


0.652 


0.610 





\Tfiv x Vp] [T ph/ x V7] 


II 


-0.153 


-0.143 






In this section we present one-loop matching coefficients 
for these operators @ 



Linear combinations of these operators are also needed for cal- 
culating the K —y 7r matrix elements of the 1 = 3/2 part of the 
electromagnetic penguin contribution to e'/e. In the context of 
staggered fermions this is explained in Ref. [24| . 



TABLE IV. Matching coefficients dj at (continued from Ta- 
ble IIII|) . Again, Greek indices are implicitly summed, with 
the condition that they are unequal, with the further con- 
straint that for the operator [A M x T up ] [Ap x T vp \ , v < p, and 



for the operators [T M „ 


x A p ] [Tfj,i, x A p 


and [Tp„ 


x Ap\ \Tpv 


xA, ; ] 


fJb < v. 










u .i 


color trace 


(c) 


(d) 


Tj 




I 


-0.026 


-0.017 





[X^y x -A^][T^t, x 


II 


-0.045 


-0.030 





[X^y X -A^^X^l, X j4,y 


I 


-0.010 


-0.007 





[T^v x -A/j,] [X^j,^ x A v 


II 


-0.003 


-0.002 





[T^v x -A^j,] [X^j,^ x Ap 


I 


0.001 


0.000 





[X^y x A.^][X^t/ x Ap 


II 


-0.002 


-0.001 





[X^y x Ap\ [Tp,i, x Ap 


I 


0.001 


0.000 





[X^ x -Ap] \Tp,i/ x Ap 


II 


-0.002 


-0.001 





[Tpv x A.p] \Tp,i/ x Ap] 


I 


-0.068 


-0.046 





[Tpv x A.p] \Tp,i/ x Ap] 


II 


-0.024 


-0.016 





[Xp^y x A.p] [Tp,i/ x Arj 


I 


0.000 


-0.000 





[Xp,y x A.p] [Xp,^ x A^ 


II 


0.003 


0.002 





L4p x S]L4p x S] 


I 


-0.003 


-0.002 





[A„ x S*]L4p x 5] 


II 


-0.022 


-0.014 





[Ap x Xp,(y][Ap, x Tp,u 


I 


-0.124 


-0.086 





[Ap, x Xp,,y][Ap, x Tpu 


II 


-0.114 


-0.084 





[Ap, x Xp,,y][Ap, x Xpp 


I 


0.022 


0.017 





[Ap X Xp,;y][Ap, X Xp,p 


II 


-0.002 


-0.002 





[Ap, x Xp,(/][Ap, x Tup 


I 


-0.004 


-0.003 





[Ap, x Xp,,y][Ap, x Tup 


II 


0.011 


0.009 





[Ap, x Xi/p] [Ap, x Tp, u 


I 


-0.004 


-0.003 





[Ap, x X,/p] [Ap, x Tp, u 


II 


0.011 


0.009 





[Ap x X^p] [Ap, x T V p] 


I 


-0.106 


-0.083 





[Ap x X^p] [Ap, x T V p] 


II 


-0.074 


-0.055 





[Ap, x T V p\ [Ap, x Turj^ 


I 


-0.017 


-0.015 





[Ap, x X^p] [Ap, x Turj 


II 


0.011 


0.009 





[PxVp][PxVp] 


I 


2.566 


2.004 


-1 


[Px Vp][Px Vp] 


II 


-0.547 


-0.503 





[P x Vp][P x V v ] 


I 


0.151 


0.141 





[PxKJPxVJ 


II 


0.326 


0.307 





[P X Ap][P X Ap] 


I 


-0.063 


-0.041 





[P x Ap][P x Ap] 


II 


-0.042 


-0.028 





[P x Ap][P x A v ] 


I 


0.012 


0.008 





[P x Ap][P x A v ] 


II 


-0.005 


-0.003 






As for the first step is to match the operators 

into the augmented continuum theory. The result is 



Cont' 



O 



^Cont' 



/7-)Cont ; 

— ^S2 

1 



U P2 ' 



- py 

2 \ Usi 



Cont' 



OCont 
PI 



/'oCont 
U T1 



1 



— 

— U S1 

1 



Cont' 



Cont' 



opt 



~ 2 { 32 



/oCont 
U P2 



o 



Cont' 
T2 



o 



Cont' 



/r-)Cont' 

— ^S2 



/oCont' 
U P2 ' 



2 I u vi 



Cont' 



U A1 



o 



Cont' 



/7-jCont' 
U S1 



OCont' , 
PI + 



/-nCont' /-/-jCont 
U V2 ~ U A2 



(41) 
(42) 
(43) 
(44) 
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where some of the operators are defined in Eqs. (|26H29I) . 
and the others are 



U S1 


= [S«(1®&)I>6][S&(1« 


>6K], 




(45) 


U S2 


= [£«(1®&)I>«][S£(1€ 


5 6)^], 




(46) 


//-jCont' 
U P1 


= [5a(75®e 5 )A][^(75 






(47) 


//~\ Coiit 
U P2 


= [S a (75® (75 






(48) 


/■/-jCont' 
U T1 


= $3[^o(7m7v®6)A] 






> 

(49) 


/■/-jCont' 
W T2 


= E^TmT^^AJ 


[^(7^ C 




• 



(50) 

At tree-level, the new operators in the augmented contin- 
uum theory match onto lattice operators in the obvious 
way: 



U S1 


tree 


/7")Lat 


= [5x P][5 


xP]/, 


(51) 




tree 


O h s f 


= [5x P][S 


x P]n , 


(52) 


/?-)Cont' 
U P1 


tree 


U P1 


= [P x P][P 


xP]i, 


(53) 


Ly P2 


tree 


l - y P2 


= [P X P][P 


X P]// , 


(54) 


/7-)Cont' 
L/ T1 


tree 




= Ep^ x 


P][T^XP]!, 


(55) 



(9 



Cont' t££ c / Lat 
T2 



= ^T2* = E^ x x P] n . 



(56) 



ll<V 



As for 0§° nt ', each of the operators 0§°f' matches at 
one loop order onto class (A) lattice operators composed 
of bilinears with taste £5, and class (B) operators having 
other tastes. Continuum mixing involves only class (A) 
operators — mixing with operators of class (B) is a lattice 
effect. We display here only results for mixing with class 
(A) operators, for several reasons. First, these are likely 
to be the subset of operators used in numerical simula- 
tions. Second, the contribution of class (B) operators is 
of next-to-next-to-leading order in staggered chiral per- 
turbation theory, using the power counting of Ref. [30j. 
This is because the contribution is suppressed both by 
a s , and by the need to have a chiral loop due to the 
mismatch between the taste of the operator and external 
states. Unlike for 0§ ont , there is no chiral enhancement 
to raise the contribution to next-to-leading order. Fi- 
nally, we do not show results for class (B) mixing for the 
sake of brevity. 

We write the one-loop matching formula as 



O 



with i 



Cont' 



E z 

- 4 and 



Lai 



(4tt) 2 ^ 



^0]^ : 



(57) 
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(47T) 



f ( 

— ^- 7 ylog(H 



jCont 



jLat 



CfImfTv 



(58) 



We stress that we are using a different basis for the con- 
tinuum and lattice operators, so that the tree-level con- 
tribution is no longer Sij. Because of this, we denote 
the finite parts as rather than Cy. The df° nt ma- 
trix comes from Ref. Q • The elements of ci^ at , which are 
calculated numerically, are new results. The TV, term 
is present only if the operator is mean-field improved. 
Numerical values for Im f are given in the caption of Ta- 
bleED 

Results for the four operators are presented in Ta- 
bles |V]|VlTTJ We show results for HYP-smeared operators 
without mean-field improvement, but include the values 
of Tij so that mean-field improvement can be easily im- 
plemented. 

TABLE V. Matching coefficients entering Eq. (JSHJ for i = 
2, i.e. for O% ont ' . The finite lattice coefficients are for the 
HYP-smeared fermion action and operators, and either (c) 
the Wilson gauge action or (d) the Symanzik gauge action 
(following the labeling used in Table U The df° nt matrix 



comes from Ref. 
digits quoted. 



Results are accurate to the number of 



Operator j 


b 2 j 


72j 


^Cont 
«2j 


dkf{c) 


d\f{d) 


T 2J 


U S1 


-1/2 


6 


-59/12 


2.70 


2.34 


-1 


U S2 


1 


-10 


+73/12 


-17.80 


-14.53 


6 


/?-)Lat 
U P1 


-1/2 


6 


-59/12 


3.64 


3.17 


-1 


/?-*Lat 
<^P2 


1 


-10 


+73/12 


5.42 


4.06 


-2 


/?-*Lat 
U T1 

/?-)Lat 
<J T2 


1/2 



-14/3 
2/3 


+29/12 
-5/4 


-2.94 
0.03 


-2.52 
0.01 


1 




TABLE VI. Matching coefficients entering Eq. ([58 
Notation as in Table IVl 



for i - 



Operator j 


b 3 j 


73i 


jCont 


jLat / \ 

d 3j (c) 


d\f (d) 


T 3j 


<^S1 


1 


8 


1/3 


-4.35 


-2.88 


2 


/^>Lat 
<^S2 


-1/2 





-5/3 


7.13 


5.38 


-3 


/?-)Lat 
Upi 


1 


8 


1/3 


-6.23 


-4.56 


2 


/?-)Lat 
{ ~ J P2 


-1/2 





-5/3 


-0.25 


-0.14 


1 


/?-)Lat 
U T1 





8/3 


-1 


0.14 


0.24 





/?-)Lat 
U T2 


1/2 


16/3 


-1/3 


-2.31 


-1.55 


1 



TABLE VII. Matching coefficients entering Eq. 
Notation as in Table IVl 



581) for i = 4. 



Operator j 


b 4j 


74j 


iCont 


d\f{c) 


d\f (d) 


Tij 


/7-\Liat 

Vsi 








-3 











U S2 


1 


-16 


+23/3 


-19.12 


-16.02 


6 


/r\Lat 
U Pl 








+3 











/r^Lat 
U P2 


-1 


16 


-23/3 


-6.93 


-5.09 


2 


/^Lat 
U V1 


-1/2 


8 


-23/6 


3.02 


2.72 


-1 


/r\Lat 
^V2 








+3/2 


0.37 


0.34 





z^Lat 
<^A1 


1/2 


-8 


+23/6 


-3.27 


-2.95 


1 


<->A2 








-3/2 


0.40 


0.37 
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TABLE VIII. Matching coefficients entering Eq. (J58J) for i = 
5. Notation as in Table IVl 



Operator j 


b 5j 


75j 


jCont 
d 5j 


d&j (c) 


d Bj (d) 




&S1 


1 


I 


1/6 


/I C7 

-4.D( 


O A O 

-o.4z 


o 
I 


^S2 





-6 


-1/2 


-2.32 


-2.45 





U P1 


-1 


-2 


-1/6 


6.55 


5.10 


-2 


/r\Lat 
U P2 





6 


1/2 


-3.32 


-2.59 





/r\Lat 
U V1 





3 


1/4 


0.16 


0.27 





/r\Lat 
U V2 


-1/2 


-1 


-1/12 


5.24 


4.04 


-2 


/r\Lat 
^Al 





-3 


-1/4 


-0.16 


-0.27 





/r^Lat 


1/2 


1 


1/12 


0.05 


0.09 






The tables show that, for coefficients with magnitudes 
larger than about 5, improvement of the gauge action 
leads to a small reduction in the size of the finite lattice 
coefficients @ Even with this improvement, we see that 
the largest lattice coefficients have a magnitude as large 
as 16, although most have a magnitude smaller than 10. 
These numbers are larger than those we found for the Bk 
operator (Table [J). On a coarse lattice with a ss 0.12 fm, 
and <7 2 /(47r) 2 ss 0.025, the largest one-loop term could 
give a 40% correction. This is not a precise statement be- 
cause we have not included the contribution from finite 
continuum terms and from the anomalous dimensions. 
Nevertheless, it is a warning to expect perturbation the- 
ory to be less convergent for the operators 02-5 than for 
Ob k - 

The situation is better if one implements mean-field 
improvement. As can be seen from the Tables, this re- 
duces all the larger coefficients, so that the largest mag- 
nitude is now « 10. Since mean-field improvement is 
relatively straightforward to implement, our results sug- 



gest that this would be worth the required investment. 

Finally, we note that, unlike for Ob k , one finds very 
large finite coefficients if one uses mean-field improved 
naive fermions (with either gauge action). The largest 
magnitudes are ~ 50, indicating a complete breakdown 
in the convergence of perturbation theory. Because of 
this, we do not include the results in the tables. 

V. CONCLUSION 

We have calculated the matching factors for four- 
fcrmion operators using various fermion and gauge ac- 
tions. Most useful are our results for the fermion action 
and operators constructed using HYP-smeared links with 
the Symanzik improved gluon action. These are needed 
for our on goin g calculation of Bk and related matrix el- 
ements [3 l3l|. For these operators, the one-loop correc- 
tions are of moderate size for the Bk operator, with the 
range of corrections being sa 10 x a s /(47r) « a s , which 
is the naively expected size. The same holds true for 
the operators induced by new physics, as long as one im- 
plements mean-field improvement. For all operators, we 
find that improving the gauge action generically leads to 
a reduction in the size of one-loop matching coefficients, 
but that the effect is relatively small. 
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